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The  Interaction  Cf  Finite  Amplitude  Deflection  And 

Stretching  Waves  In  Elastic  Membranes  And  Strings. 

By  D.  F.  Parker 

Department  of  Theoretical  Mechanics^ 

The  University,  Nottingham. 

an  d 

E.  Varley 

Center  for  the  Application  of  Mathematics^ 

Lehigh  University. 

SUMMARY. 

Disturbances  produced  by  the  motion  of  a  driver  which  is 
rigidly  bonded  to  the  edge  of  a  plate  are  used  to  motivate  parameter 
expansion  techniques  which,  when  applied  to  the  equations  of  finite 
elasticity,  generate  approximating  equations  which  describe  low 
frequency  deflection  and  stretching  waves  travelling  along  stretched 
elastic  plates  and  rods  in  the  limit  when  bending  forces  are 
negligible  compared  with  membrane  forces.  The  structure  of  the 
boundary  layer  at  the  driven  edge  and  shock  layers  ,  where  the  low 
frequency  or  filament  approximations  are  locally  invalid,  are  also 
discussed. 

The  low  frequency  equations  are  used  to  discuss  the  interaction 
between  progressing  finite  amplitude  deflection  and  stretch  waves 
in  the  limit  when  the  stretch  rate  is  small  compared  with  the  angular 
speed  of  the  plate.  The  disturbance  is  locally  that  of  a  pure 
deflection  simple  wave  whose  amplitude  and  frequency  are  modulated 
by  slow  variations  in  the  stretch.  As  the  stretch  increases  the 
frequency  increases  while  the  amplitude  decreases.  The  stretch  vave 
is  also  modified  by  deflection  of  the  plate:  the  speeds  of  wavelets 
carrying  constant  values  of  stretch  are  always  less  than  their 
values  in  the  pure  stretch  simple  ware. 
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1 .  Introduct i on. 

This  paper  is  in  two  parts.  In  the  first  part  we  describe 
an  expansion  procedure  which,  when  applied  to  the  equations  of 
finite  elasticity,  generates  equations  which  govern  low  frequency, 
finite  amplitude  ,  deflection  and  stretching  disturbances  in 
elastic  plates  in  the  membrane  limit  when,  to  a  first  approximation, 
bending  forces  are  negligible  compared  with  membrane  forces.  The 
first  term  in  the  expansion  giveathe  classical  membrane,  or  flexible 
string  approximation,  (see  Cristescu  (5)  and  (6),  and  Craggs  (7)). 

In  the  second  part,  these  low  frequency  equations  are  used  to 
discuss  the  interaction  between  finite  amplitude  deflection  pulses 
and  stretching  waves.  Detailed  analysis  is  given  for  plates  and 
is  directly  applicable  (see  Parker  [ 10  ]  )  to  stretched  rods  in  the 
'flexible  string'  limit  when  bending  moments  and  torques  are 
negligible . 

A  region  is  a  low  frequency  region  for  a  deformation  if: 

(I)  the  strain  rate  and  angular  velocity  of  the  plate  at  all 
particles  are  (negligibly)  small  compared  with  the  frequency  defined 
by  the  local  sound  speed  C  and  plate  thickness  d; 

(II)  the  plate  curvature  and  stretch  gradients  at  all  particles 
are  small  compared  with  d 

In  section  2  the  expansion  procedures  are  motivated  by  considering 
the  disturbances  which  are  produced  when  one  edge  of  the  plate  is 
rigidly  bonded  to  a  driver  which  deflects  and  stretches  the  plate. 

It  is  supposed  that  the  angular  velocity  and  stretch  rate  at  the 
edge  of  the  plate  are  proportional  to  w^d/C  and  Wgd/C:  for  time 
periodic  motions  these  can  be  taken  as  the  frequencies  defined  by 
the  variations  in  direction  and  magnitude  of  the  applied  traction. 

The  dependent  and  independent  variables  are  then  non- di mens i onal i b ed 
in  such  a  way  that,  in  terms  of  these  variables,  the  governing 
equations  have  formal  regular  asymptotic  solutions  which  are  valid 
as  w=max{u>  ,w  }-*0.  These  satisfy  the  traction  free  conditions  on 
the  lateral  surfaces  of  the  plate,  and  describe  deformations 
satisfying  conditions  I  and  II.  The  lowest  order  approximation 
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agrees  with  the  classical  filament,  or  long  wave,  approximation  and 
is  analogous  to  the  long  gravity  wave  approximation  used  in  hydraulics. 

The  filament  theory  cannot  satisfy  the  bonding  condition  imposed 
at  the  edge  of  the  plate  since  this  is  inconsistent  with  the  Poisson 
contraction  which  occurs  in  the  low  frequency  region.  A  boundary 
layer  (called  the  first  diameter  by  experimentalists  (see  Bell  [l]) 
occurs  in  which  conditions  II  are  violated.  In  section  3.2  we 
introduce  boundary  layer  variables  and  show  how  the  low  frequency 
expansion  must  be  modified.  To  a  first  order  inertia  effects  are 
negligible  and  the  applied  tractions  and  particle  velocities  are 
continuous  across  the  layer. 

If  the  plate  is  not  deflected  the  low  frequency  approximation 
neglects  the  lateral  inertia  of  the  plate  compared  with  its 
longitudinal  inertia.  In  section  3.1  we  use  the  well  known  pure 
stretching  simple  wave  solutions  to  show  how  the  strain  rate  and 
lateral  particle  speeds  can  grow  until  the  low  frequency  approximation 
becomes  invalid.  A  'shock  layer'  forms  in  which  both  sets  of 
conditions  I  and  II  are  violated:  in  this  layer  the  lateral  and 
longitudinal  inertias  of  the  plate  are  comparable. 

In  section  3.3  we  exhibit  solutions  of  the  filament  equations 
which  describe  pure  deflection  simple  waves;  these  deflect  but  do 
not  stretch  the  plate.  Such  progressing,  non-distorting  waves 
cannot  form  shocks.  The  second  half  of  the  paper,  section  ^  , 
emphasises  the  role  played  by  these  waves  when  the  stretch  rate  is 
small  compared  with  the  angular  speed  of  the  plate.  Such  deformations 

occur  when  e  =  w  / oj  <<  1. 

S  D 

The  interaction  between  a  deflection  and  stretch  wave  is 
illustrated  in  section  U.l,  where  discontinuity  analysis  is  used  to 
determine  the  variation  in  the  angular  velocity  of  the  plate  at  a 
front  which  separates  the  interaction  region  from  a  pure  stretch 
precursor  wave.  If  the  stretching  wave  extends  (contracts)  the 
plate  as  it  passes  the  front,  the  angular  speed  increases  (decreases). 


In  particular,  when  the  plate  is  completely  unloaded  the  angular 
speed  tends  to  zero,  although  the  curvature  of  the  plate  increases 
without  bound  until  conditions  II  are  violated  and  bending  forces 
become  important. 

In  sections  h.2  and  k . 3 ^ mot i vat e d  by  the  results  of  section 
U.  1  and  previous  work  on  pulse  propagation  by  Varley  and  Cumberbatch 
(13)  and  (lU),  and  Varley  and  Rogers  (15),  we  present  expansion 
techniques,  which  are  valid  as£-vO,  to  describe  conditions  in  an 
interaction  pulse  where  the  plate  is  being  deflected  much  moie 
rapidly  than  it  is  being  stretched.  No  restriction  is  placed  on 
the  amplitude  of  the  disturbance,  and  the  'shapes'  of  the  stretch 
and  deflection  profiles  are  quite  arbitrary.  Even  though  the 
expansions  only  describe  conditions  at  a  particle  over  a  limited 
time,  no  restriction  is  placed  on  the  distance  the  pulse  propagates 
and  the  modulation  of  high  frequency  deflection  pulses  by  slowly 
varying  stretching  waves  is  described  in  detail.  The  deformation, 
at  any  time,  is  locally  that  of  a  deflection  simple  wave  with 
amplitude  and  frequency  modulated  by  the  slow  variations  in  stretch. 
Its  amplitude  increases  as  the  local  stretch  decreases  and  becomes 
unbounded  when  the  tension  is  completely  released.  The  stretching 
wave  is  itself  affected  by  the  mean  values  of  the  deflecting 
disturbance.  Wavelets  carrying  constant  values  of  stretch  are 
always  slowed  relative  to  their  speeds  in  the  pure  stretch  wave. 

Finally,  in  section  5,  ve  discuss  expansion  techniques  for 
deformations  in  'shock  layers'  in  which  both  sets  of  conditions 
I  and  II  are  violated.  In  general  the  speeds  of  these  layers  are 
not  those  of  thermodynamic  shocks:  dissipation  is  negligible  and 

the  deformation  is  isentropic.  The  effect  of  a  shock  layer  on  the 
low  frequency  deformation  can  be  obtained  by  treating  the  layer  as 
a  shock  discontinuity  across  which  conservation  laws  can  be  applied 
to  derive  jump  conditions.  These  laws  are  identical  to  those  defined 
by  writing  the  filament  equations  in  appropriate  divergence  form. 

If  a  shock  layer  separates  two  low  frequency  regions  with  identical 
deformations  at  the  edge  of  the  layer,  conditions  in  the  layer  are 
identical  to  those  in  a  solitary  wave. 
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2 .  General  formulation. 

Our  aim  is  to  construct  a  model  which  describes  the  main 
features  of  the  deformations  produced  in  stretched  elastic  plates 
and  rods  by  'low  frequency'  deflecting  and  stretching  waves  which 
are  controlled  by  the  presence  of  traction  free  boundaries. 

Detailed  analysis  is  presented  for  plates,  and  the  corresponding 
results  for  rods  quoted. 

The  deformation  is  described  in  the  usual  way  by  explicit 
relations 

x  =  x( X , t )  (2.1 ) 

between  the  cartesian  co-ordinates  f  of  a  particle  at  time  t  and 
its  co-ordinates  X  in  some  reference  configuration  R  where  the 
plate  occupies  the  region  X^O,  -d£X2£d,  and  -®<X3<®,  and  is  in 
equilibrium  under  uniform  hydrostatic  pressure*,  uniform  temperature, 
and  uniform  density  0.  We  consider  adiabatic  deformations,  which 
occur  when  the  power  generated  by  the  heat  flux  or  other  energy 
sources  is  negligible  compared  with  the  rate  at  which  stresses  do 
work,  and  restrict  our  attention  to  elastic  materials  which  are 
isotropic  and  homogeneous  with  respect  to  their  states  in  R.  For 
such  a  material  the  internal  energy  is  a  function  of  the  entropy 
S  and  any  three  i  n  depe  nden  t  invar  i  ant  s  of  the  Cauchy-Green  strain 
tensor 

G  =  pTp ,  (2.2) 

where  p,  the  deformation  gradient  tensor,  has  components 

ax. 

P  L  j  =  1/3Xj-  (2.3) 

We  con  r  plane-strain  deformations  for  which 

Xj  Uj  ,X2  ,t  ),  x2  =  x2(Xj,X2,t)  and  x3  =  X3,  (2.U) 

*  T) T 3  aken  as  the  ambient  pressure.  In  what  follows  all 

st-»3s#s  are  measured  relative  to  this  hydrostatic  stress. 
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which  for  t>0  are  isentropic  with  S  =  Sq.  As  the  three  invariants 
of  G  we  take 


I  = 


2  2  2 

2 ( t  r . G- 1  )  ,  J  =  (det.G)5,  and  K  =  2I  +  J  -|[(tr.G)  -tr.g‘):  (2. 


5) 


in  plane  strain  deformations  K=0. 

The  elastic  material  is  specified  by  its  internal  energy 
oE(l,J,K:S)  ,  where  o  is  a  characteristic  constant  of  the  material 
with  dimensions  of  stress,  (e.g.  Young's  modulus  at  ambient  strain), 
and  E  is  non-dimensional.  If  we  define 


W(I,J,S0 )  =  E ( I , J , 0 : S  o ) 


(2.6) 


then,  for  plane  strain  isentropic  deformations,  (see  Green  and 
Adkins  (8),  p.lOU),  the  two  dimensional  Piola-Ki rchof f  stress 
tensor,  n on - di mens i on al i s ed  by  o,  is 


T  =  (TiJ) 


.fp" 

VP2  1  P2  2  J  y"1 


22 


1  2 


where  henceforth  the  subscripts  take  the  values  1,2. 

It  is  convenient  to  work  with  a  system  of  first  order  equations 
for  T,  p  and  the  particle  velocities 

3  x .  . 

u.=  13t.  (2.8) 


Consequently,  to  equation  (2.7),  and  the  Euler  equations 


3T 


iJ/9X  =  C'2  Sui/3,, 
J 


(2.9) 
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where  C  =  /o  /  P  is  a  characteristic  speed  defined  by  the  material 
we  add  the  compatibility  conditions 
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which.,  together  with  suitable  initial  conditions,  imply  (2.3) 
and  (2.8).  During  the  deformation  the  surfaces  of  the  plate 
X2*+d  are  traction  free  so  that  (2.7),  (2.9)  and  (2.10)  are  to 
be  solved  subject  to  the  conditions 


T22  =  T12  =  0,  on  X2  =  +d. 


(2.11) 


The  disturbance  in  the  plate  is  produced  by  the  forced 
motion  of  its  edge  X^O.  For  definiteness  we  take  this  edge 
to  be  rigidly  bonded  to  a  moving  plate  so  that 

p22-1  =  PI2  =  0.  on  X,  =  0.  (2.12) 


This  driver  transmits  a  total  force  per  unit  breadth  to  the 
plate  with  varying  magnitude  F(t)  an-d  di  re  ct  i  on  e(  t  ) .  The 
variation  of  F(t)  stretches  the  plate,  and  the  variation  in 
q(t)  deflects  it.  We  suppose  that  the  time  variations  of  ?  and 
6  introduce  two  characteristic  frequencies  an<i  wD  so  that 

F'(t)  =  2dousF'  (  isT)  ,  and  0 '  (  t  )  =  “D  e'  (  “DT)|  (2.13) 

where  F  and  0  are  bounded  functions  of  their  arguments  with 
piecewise  continuous  bounded  derivatives.  In  the  special  case 
of  time  periodic  variations  in  applied  force  and  can  be 
taken  as  the  frequencies  defined  by  the  variations  of  F  and  6. 

In  order  to  make  the  statement  'low  frequency  wave'  precise 
we  work  with  non-dimensional  independent  variables 

t  =  ill x  ,  X  =  “X  /C,  and  Y  =  X2/d  (2.1*0 
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where  u>=  max .  {  fig  ,  u>D }  .  As  non-dimensional  dependent  variables 
we  take  T,  p  and 


(v, 


|V2} 


-  1  / 

C  ( u. 


u2  ) 


(2.15) 


In  terms  of  these  variables  conditions  (2.7)  are  unchanged; 
(2.9)  are  replaced  by 


3T 


1  2  / 


3  Y 


3  V 

t»)C 


1 


(2.16) 


an  d 


3T 


22/ 


3  Y  =  u) 


r> 


/ 

3t  - 


3T 


2  l/ 


3X), 


(2.17) 


whe  re 


Q)  —  U)  d  /  C  y 


(2.18) 


and  (2.10)  are  replaced  by 


3p 


11/ 


3  v 


3 1  - 


1/ 


3  X  = 


3  P  2  1  / 


3  V 


3 1 


2  / 


3X  =  0 


(2.19) 


an  d 


*'■>/„.  ■  »3P22/3X.  (2.20) 

The  parameter  w  occurring  in  ( 2 . 16 ) - ( 2 . 20  )  in  the  ratio  of 
the  maximum  frequency  of  the  traction  imposed  at  the  edge  of  the 
plate  to  the  characteristic  frequency  C/d  defined  by  the  material 
and  the  stress  free  boundaries.  Alternatively,  it  is  the  ratio  of 
plate  thickness  to  a  typical  wavelength  of  the  disturbance. 
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Here  ve  are  interested  in  low  frequency,  long  wave,  disturbances 
for  which  u><<1.  More  precisely,  a  region  is  a  low  frequency  region 
if  as  w+O  the  functions  T(X,Y,t:w),  g(X,Y,t;u>)  and  v(X,Y,t:<*>)  have 
limiting  values,  satisfying  the  stress  free  conditions  (2.11).  These 
functions  have  bounded  derivatives  which  need  nbt  be  identically  zero. 

- g - - - - 

If  f  denotes  the  limiting  value  of  any  variable  f  as  w-*o  then 
(2.11),  (2.16)  and  (2. IT)  imply  that 

T12  h  T22  e  0,  (2.21) 

so  that,  by  (2.7),  the  components  of  p  are  implicitly  related  by 
the  conditions 


o  o  o  o  .  . 

pp  +PP  =0  (2.22) 

11  12  2221 


and 


0  3W/„t  a  0 

p  3 1  +  p 

22  11 


3  W/ 


3  J  =  0. 


(2.23) 


(2.19),  (2.20),  (2.22)  and  (2.23)  then  imply  that  3p/3Y  =  0  so 
that 

p  =  p (X , t ) .  (2.2k) 

Conditions  (2.22)  and  (2.23)  imply  that  in  the  limit  of  low  frequency 

disturbances  the  families  of  fibres  X  =  constant,  Y  =  constant  which 

* 

are  mutually  orthogonal  in  H^remain  orthogonal  and  in  directions  of 

principal  stretch.  Moreover,  like  the  free  surfaces  Y  =  ±1,  any 

surface  Y  =  constant  is  stress  free.  If  (A,p)  denote  the  limiting 

o  o 

principal  stretches  and  (T  T  )  the  principal  Pi ol a-K i r ch o f f  stresses 
then,  in  terms  of 

G  ( A  ,A  :S  )  =  W  (  2  (  A  2  ♦  A2)  ,A  A  :S  )  , 

12  0  1  2  12  0 


(2.25) 


[10] 


A  and  p  are  implicitly  related  by  the  condition 


o 

T 

2 


U,w:S0  ) 


0, 


(2.26) 


while 


O  3  Q 

T1  =  TT  »H : S0  )  *  =  T(l )  say* 
3A1 


(2.27) 


where  the  dependence  of  T  on  S0  is  not  stated.  Relation  (2.27)» 
with  u(A)  given  by  (2.26),  is  exactly  the  relation  between  force 
per  unit  undeformed  area  (engineering  stress)  and  extension  for 
simple  static  stretching  of  the  plate  with  ambient  temperature 
©adjusted  so  that 

©=  -  °T?(  I»J»0;S0).  (2.28) 


3 . 1  Asymptotic  expansion  for  the  low  frequency  region. 


In  the  low  frequency 

region 

p ,  T  and  v  admi t 

asympt  ot i c 

expans  ion 

p  =  |(x,t) 

N 

r  n 

♦  l  <*» 

p(X,Y,t), 

n=5 1 

T  =  T  ( X  ,  t ) 

N 

♦  l  *n 

?(X,Y,t), 

n=  1 

and 

v  =  v  ( X  ,  t ) 

N 

c  n 

+  1  u 

v  ( X  ,  Y  ,  t )  . 

n=l 

When 

(3.1) 

is  inserted,  (2 

.  l6  )  -  (  2 

. 20 )  imply  that 

’*>'«  -  3?"' 

o 
3  v 

at  = 

2  /  ^2  1/ 

'  ax  -  W3t  = 

and 

'  at 

.  ”»/„ 

c 

3V2/ 

at 

formal 


(3.1) 


(3.2) 


(3.3) 


Since  the  right-hand  sides  of  (3.3)  are  independent  of  Y,  and  since 


[12] 


and 


3$ 


1/ 


3$ 


at  - 


11/ 


o 

3  v 


3  X  = 


2  / 


3 1  - 


2  1 


3  X  =  0. 


(3.5) 


Equations  (3.2)  and  (3.5),  together  with  (2.7),  (2.22)  and  (2.23), 
provide  a  complete  set  of  equations  for  the  zeroth  order  variables. 
They  are  the  low  frequency  membrane  equations  which  are  usually 
derived  by  postulating  (2.21)  and  assuming  p  is  independent  of  Y, 
(see  Craggs  (7)).  These  conditions  are  more  conveniently  written  if 
we  introduce  as  basic  dependent  variables  the  angle  between  the 
fibres  Y  =  constant  and  the  x^axis,  and  (v,w)  the  components  of 
particle  velocity  along  and  normal  to  these  fibres.  In  terms  of 
these  variables,  the  stretch  X,  and  the  tension  T(x), 


3*/. 


T  '  ( X )  '  3X  =  at  -  w  ^ '  3 1 , 


Sx/St  =3V/3X  - 


(  3.  7) 


an  d 


T  ( X  )34'/3X  =  3w/3t  +  v3^/3t, 

(3.8) 

X3^3t  =  3  w/3  X  +  v3,P/3X. 

Equations  (3.7)  and  (3.8)  are  the  equations  governing  the  plane 
vibrations  of  an  extensible  string  with  T(x)  representing  the  tension 
in  the  string.  They  also  govern  the  one  dimensional  motion  produced 
when  an  infinite  half  -space  of  incompress ible  elastic  material  is 
sheared  at  X=0  by  a  force  which  varies  in  strength  and  direction  with 
time  (see  Chu  (2),  and  Collins  (3)):  then  X  represents  the  shear  strain^ 

T(x)  the  shear  stress,  \Ji  the  direction  of  shear,  and  (v,w)  the  components 
of  velocity  parallel  and  normal  to  the  shear  direction. 
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We  restrict  attention  to  situations  when  the  plate  is  in  tension 
(T>0)';  when  T<0  bending  forces  are  important.  When  T>0  equations 
(3.7)  and  (3.8)  are  totally  hyperbolic  with  characteristic  Lagrangian 
speeds 

1 

dX/dt  =  ±  [  T '  ( A  )  ]2  ,  =  ±  Cg  (  A  )  say,  (3.9) 


and 


1 

dX/dt  =  ±  [T(a)/a]4,  =  ±  Cd(a)  say.  (3.10) 

If  equations  (3.7)  and  (3.8)  are  formally  linearized  about  A  =  AQ> 
and  v=w=  iJ/  =  0  the  pair  (3.7),  which  describes  the  stretching  of  the 
plate,  uncouple  from  the  pair  (3.8),  which  describes  its  deflection. 
(A,v)  and  (ij>,w)  then  satisfy  the  linear  wave  equations  with  constant 
Lagrangian  speeds  Cg(AQ)  and  0^(1^).  In  particular,  the  linear 
theory  predicts  that  any  disturbance  entering  a  uniform  region  which 
Is  prestrained  by  a  tension  T ( A  q )  is  generated  by  two  progressing, 
undistorted,  noninteracting  waves:  a  longitudinal  wave  composed  of 

wavelets,  each  moving  with  the  same  Lagrangian  speed  Cg(AQ),  at  which 
(a,v)  are  constant,  and  a  transverse  wave  with  wavelets  moving  with 
the  same  Lagrangian  speed  CD(AQ)  at  which  (\p,v)  are  constant.  In 
terms  of  the  driving  conditions  (2.13) 

T  =  T(a0)  +  FUgft  -  Cg 1 X  ]  )  ,  and  *  =  0  ( uD  [  t  -  C^X]),  (3.11) 


whe  re 

<»>s  =  ig/oi,  and  =  u>D/u>.  (3.12) 

The  only  known  exact  solutions  to  the  non-linear  equations  (3.7) 
and  (3.8)  which  describe  disturbances  moving  into  an  undisturbed 

region  are  the  simple  wave  solutions.  These  are  of  two  dinstinct 


types.  Stretching  (longitudinal)  simple  waves  for  which  (A,v)  vary 
while  are  constant,  so  the  plate  is  not  de f 1 ecte d^ an d  deflection 

(flexural)  simple  waves  for  which  (  ip ,  w )  and  v  vary  while  A  is  constant, 
so  the  plate  is  not  stretched.  We  might  suspect  that  these  also 
provide  a  good  'local'  description  of  deformations  for  which  ios/u>D>>l 
and  a)g/u)p<<l.  In  sections  3.1  and  3-3  we  describe  these  simple  waves 
and  in  section  U  compare  these  with  asymptotic  solutions  to  (3.7)  and 
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3  •  2  Stretching  simple  wave. 

The  stretching  simple  wave  is  governed  by  the  non-linear 
equat ions 


(X  )aX/3X 


-  9  v/  =  3  X  / 


9  t 


at 


3\ 


ax 


=  0 


(3.13) 


which  are  discussed  in  detail  by  Courant  and  Friedrichs  (4). 

For  such  waves  (X,v)  are  constant  on  wavelets  which  travel  with 
constant  but  distinct  Lagrangian  speeds.  If  fj(X,t)  =  constant 
is  the  wavelet  which  left  X  =  0  at  t  =  6 ( >0  )  then 

X  =  A(g)  (  3 . 1*0 


is  arbitrary; 


X  =  Ug(0)(t-0)  ,  (3.15) 

whe  re 

Us(6)  =  Cg ( A  )  ,  (3.16) 

an  d 

v  =  v( X  )  =  _  / cg(r)dr,  =  V ( 6 )  say.  (3.17) 

The  roles  of  these  simple  waves  are  well  understood.  Any 
disturbance  produced  by  a  driver  motion  for  which  0'  e  0  is 
generated  by  a  stretching  simple  wave  with  A(B)  determined  from 
the  condition 


T  (  X  )  =  T(  '  )  +  F  (  u>  B  )  . 

0 


A 


(3.18) 
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If  O'  -  0,  but  only  strains  satisfying 

T'(a)  >  T/a,  so  Cs(a)  >  Cd(a),  (3.19) 

are  produce  then  ,  at  any  time,  the  disturbed  region  can  be 
separated  into  two  distinct  regions!  a  region  neighbouring 
X  =  0  in  which  <1* ,  X  %  v  and  w  may  all  vary,  and  preceding  this 
a  simple  wave  region  in  which  (4%w)  are  constant  while  (A,v) 
vary  according  to  (3-1^*)— (3.  IT)  for  some  a(b)> 

In  the  simple  wave  region  the  displacement  field  is  given 
by 

£ 

ux  /d  =  A(B)X  +  V( 6)t  +  / rA ' ( r  )  U  ( r ) dr  ,  and  x  /d  =  n(A)  Y  (3.20) 
1  1 

where  p(X,t)  is  given  implicitly  by  (3.15).  In  particular  the 
equations  of  the  free  surfaces  Y  =  ±  1  are  given  by 

x2/d  =  ±  u(A).  (3.21) 

Conditions  (3.20)  imply  that  to  a  first  order  in  the  low 
frequency  region 


v 


1 


P'U) 

vTaT" 


3  X  /  v 

so  that  in  regions  where  9t  is  bounded  2/vj=  0  ( u )  as  u>-*0. 

However,  as  is  well  known,  (3.11*)  and  (3.15)  predict  that  either 

all  compressive  or  all  expansive  wavelets  begin  to  coalesce  and 
3  A  / 

dt  may  become  unbounded.  In  section  5  we  show  that  wavelets 

are  in  fact  prevented  from  coalescing  by  the  formation  of  a 

'shock  layer'  in  which  the  low  frequency  expansion  (3.1)  is  not 

3  A  /  .  . 

valid.  In  this  layer  3t  =  0(1;  as  w+O  and  the  lateral  and 

longitudinal  inertias  of  the  plate  are  of  the  same  order. 


[  IT] 


In  section  U  we  discuss  how  a  'rapid'  long  wavelength 
deflection  disturbance  modifies  a  stretching  simple  wave  of 
much  lower  frequency.  Such  disturbances  are  produced  at  X=0 
when  F ( t )  varies  slowly  compared  with  0 ( t  )  .  The  front  sep¬ 
arating  this  disturbance  and  the  pure  st-etch  simple  wave  is 
necessarily  a  characteristic  wavelet  moving  with  Lagrangian 
speed 


CD ( A ) ,  =  U  D ( 8 )  say.  (3.23) 

If  o<.(  X  ,  t  )  =  constant  is  any  such  wavelet  moving  in  a  direction 
of  increasing  X  then  (3.1*0,  (3.15)  and  (3.23)  imply  that  5: 
can  be  chosen  so  that 

X  =  xo(B)-oxi (b). 


and 


(3-2*0 


t  =  t  (pj-fftj (e). 


vhe  re 


x  -u  ( t  -e )  =  x  -  u  t  =  o  , 

0  S  0  1  SI  > 


(3.25) 


and  ( t  » t  j  )  satisfy  the  ordinary  differential  equations 


'W4; ♦  usto-(eus1'  *  (us-uD)t; +  us‘. =  °-  (3-26) 


subject  to 

t  (0)  =0,  t  (0)  =  1. 


(3.27) 


[18] 


3 .  3  First  diameter  effect 

In  addition  to  the  shock  layer  region,  where  the  lateral  and 
longitudinal  inertias  of  the  plate  are  of  the  same  order,  the  low 
frequency  expansion  is,  in  general,  invalid  in  some  vicinity  of 
the  driven  edge  X=0.  There  the  imposed  conditions  are  not,  in 
general,  compatible  with  (2.22)  and  (2.23).  Fcr  example  the  pure 
stretch  simple  wave  solutions,  given  by  (  3  •  12  )  -  (  3.  18  )  ,  satisfy 
the  boundary  conditions  (2.13)  with  6'  =  0,  but  produce  varying 

strain  fields  which  do  not  satisfy  the  bonding  condition  (2.12) 
which  requires  that  P22=l  at  X  =  Q.  In  fact  a  boundary  layer  in 
which  X=Q(u))  occurs  vhere  the  stress  and  strain  gradients  are 
larger  than  those  in  the  low  frequency  region  by  a  factor  0 ( w  ). 
In  this  boundary  layer,  which  is  called  the  first  diameter  by 
experimentalists  (see  Bell  (l))the  deformation  depends  upon  the 
precise  details  of  the  loading  distribution  on  X=0. 


In  the  boundary  layer  region  the  variations  of  the  dependent 
variables  in  the  Xj  and  X2  directions  balance.  If  we  introduce 
the  boundary  layer  variable 


-  _  1 
X  =  u  X 


X  /d 
1 


(3.28) 


equations  (2.16),  (2.17)  and  (2.20)  become 


3 T  ,  3T  ,  3  v  . 

M/3X+  1  2  ^  3  Y  -  W  1  ^  3 1  ’ 


3  T 


2  1  ^3X 


3  T 


2  2  / 


3  Y 


3  v  , 
w  2  / 


3  t 


( 3. 29  ) 


an  d 


3  P 


11/ 


_  8P12  /  _  =  3P2  j  /  _  3P 

a  Y  '  3X  1  '  3Y 


22/ax  * 


0, 
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where  T  =  T(p)  is  given  by  (2.7)-  Equations  (2.19)  yield 


3vl/.-  «  ,.3Pll/,+  .  and  3V2/9x  =  W3P21/--  (3*30) 


3  X  w 


3 1 


at " 


To  the  zeroth  order  inertia  effects  are  negligible  and  the 
only  time  dependence  is  introduced  by  the  boundary  conditions 
at  X  =  0  and  the  usual  matching  conditions  (see  Van  Dyke  (13)) 
which  require  that  as  X-*«D 

N 

v  $■  (  0  ,  t  )  +  l  uin  v  (  0  ,  Y  ,  t  ) 
n=  1 

and  (3.31) 

N 

p  -  p  (  0  ,  t  )  +  l  ajn  p(0,Y,t). 


Equations  ( 3 . 30  )  and  (3.31)  imply  that  to  the  zeroth  order 
v=v(0,t)  throughout  the  layer,  while  (3.29)  imply  that  the  stresses 
satisfy  the  equilibrium  equations 

9T  .  3  T  .  3T  ,  .  3  T  .  _  n  ( ?  \ 

u/ax+  1 2  7  a  y  2 1 7ax  22/aY"°-  (3.32) 

The  boundary  layer  deforms  like  a  semi -in  finite  plate  rigidly 
bonded  at  X  =  0  under  a  load  2doT(i)  [cosij;,  sin^]  applied  at 
X  =  °°.  This  static  elastic  deformation  has  not  been  quantitatively 
determined.  However,  because  of  the  divergence  form  of  (3.29)  the 
total  load  at  time  t  is  transmitted  across  this  layer.  Consequently, 
this  layer  does  not  affect  either  the  velocity  or  load  transmitted 
to  the  low  frequency  region,  and  conditions  (2.13),  or  equivalent 
conditions  on  the  velocity  field,  can  be  specified  as  boundary  data 
on  X  =  0  for  equations  (3.7)  and  (3.8). 


3.*<  Deflection  simple  wave. 

The  deflection  simple  wave  has  not,  to  our  knowledge,  been 
discussed  previously.  X  does  not  vary,  so  the  membrane  is  not 
stretched,  while  i|/ ,  w,  and  v  are  constant  on  wavelets  3f(X,t)  = 
constant  which  move  with  the  same  Lagrangiar  speed  C^(A).  Since 
U'  U'  5  0,  (3.210  -  (3.27)  integrate  to  give 

o  U 


t  = 


(3.33) 


or,  eliminating^ 


V 


-  X 


-  <VV“' 


The  deflection  angle 


( 3. 3k) 


is  arbitrary,  while 


ip  =  lff(  a) 


(  3. 35  ) 


v  +  (XT)  =  q  cos  Y  and 
w  =  -qs  in4>-Y 


(  3. 36) 


where  q  and  if  are  constant  Riemann  invariants.  Conditions  (3.35) 
imply  that  to  a  first  order  the  velocity  field  is  given  by 

Vj  =  -  (  A  T )  2  cos  tjj  +  qcosy  ,  and  v?  =  -(ATj^siniJi  +  qsinY,  (3.37) 

which  produces  a  displacement  field 

i.  -  x2  / 

ojX  j  /d  =  tqcosY  -  (  A  T  )*  /“  cost  (  r  )  dr  ,  and  U)(  /  d  -  y(A)Y)  =  tqsiny 


-  (  AT)  I /  i  n  y ( r ) dr . 


( 3. 38) 


If  the  simple  wave  is  moving  into  s  uniform  region  where  ^=0 
then 


q 


and  |  =0 


so  that  conditions  (3.36)  become 


1 

Vj  =  (^Tj^d-cosiJi),  and  v^  =  -(AT)^sin^. 


(  3.  39  ) 


(  3.  U0  ) 


De  fleet! on  si mpl e  waves  can  occur  adjacent  to  a  uni  form 
region  which  is  prestrained  so  that  C  >C  .  They  completely  define 
the  disturbance  in  a  membrane,  or  string,  which  is  produced  by 
changing  the  direction  but  not  the  magnitude  of  the  applied 
traction.  Then  F'  E  0  and 

U<5/ 

'P(ar)  =  e[  (  b/uD  -  l  )w] 

We  now  show  how  the  statements  (3.33)  -  (3.^0),  which  are  exact 
when  A  is  constant,  may  be  modified  to  provide  a  good  description 
of  disturbances  in  which  A,  q  and  y  are  ’slowly'  varying  compared 
with  ip.  Suon  distrubances  are  produced  when  u>c/u>  <<1. 

LJ  D 
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k . 1  Modulation  of  finite  amplitude  deflection  pulses  b 


stret  chin/ 


slow 


The  deformations  described  in  sections  3.2  and  3.^  correspond 
to  exact  solutions  of  equations  (3.7)  and  (3.8).  They  are  generated 
by  special  motions  of  the  edge  X=Q .  Although  the  pure  stretch  waves 
are  produced  quite  naturally  by  simple  rectilinear  motions  of  the 
drive: ,  the  pure  deflection  waves  are  produced  only  by  driving  motions 
which  rotate  but  do  not  stretch  the  membrane  at  X=0.  These  edge  motions 
produce  no  variation  of  stretch  at  other  values  of  X  and  are  somewhat 
artificial.  If,  however,  we  consider  those  driver  motions  for  which 
F  ( t  )  varies  slowly  compared  with  0(t)  (  S/w£j<<l)  an(i  which  produce 

progressing  disturbances  such  that  the  stretch  X,  at  some  station 
X=  X  o  at  time  t  =  t0  is  varying  'slowly'  compared  with  \p ,  it  is  not 
unreasonable  to  expect  that  in  some  interval  around  (X0*t0)  the 
disturbance  is  'approximately'  described  by  a  simple  wave  solution 
(3.35)  and  (3.36)  with  appropriately  ch<-  en  Riemann  invariants  X,  q 
and  y.  The  problem  then  is  to  construct  a  procedure  for  'enveloping' 
these  local  simple  waves  to  obtain  a  global  statement  for  conditions 
in  a  progressing  wave  in  which  X,  q  and  y  vary  slowly  compared  with  . 

Of  course  we  cannot  expect  any  one  simple  procedure  to  apply  for  all 
F  ( t )  and  6(t)  at  all  (X,t).  To  illustrate  these  procedures  we 
investigate  conditions  in  an  interaction  pulse  which  is  preceded  by 
an  arbitrary  pure  stretch  simple  wave  described  in  section  3.2. 

We  construct  formal  asymptotic  solutions  to  (3.7)  and  (3.8)  as 

w  S  / 

c  «=  -*•  0  which  describe  how  a  pure  deflection  wave  is  modulated  by 

slow  variations  of  the  stretch  X.  Even  though  the  solutions  are  valid 
at  a  given  X  only  for  a  finite  time  interval  after  the  passage  of  the 
interaction  front  a=0  of  (3.2^),  no  restriction  is  imposed  on  the 


distance  this  pulse  propagates.  Conditions  in  the  pulse  are  related 


to  6(t)  and  to  conditions  in  the  pure  stretch  wave  which  precedes  it: 


at  X=0  the  pulse  lasts  only  while  X  changes  by  0(e)  although  0  can 


change  by  any  amount. 
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The  interaction  is  best  described  in  terms  of  new  dependent  and 
independent  variables.  As  dependent  slow  'Riemann'  variables  we  take 
q  ,  y  an  d 

h 

a  =  A  C  D  =  (AT) 
and  as  our  fast  variable 


$  =  <P  -  Y  • 


(U.l) 


These  are  related  to  v  and  w  as  in  (3.36),  where  the  slow  variables 
are  Riemann  invariants.  As  independent  variables  we  take  character¬ 
istic  parameters  (a,n)  where  the  a(X,t)  =  constant  and  n(X,t)  = 
constant  wavelets  propagate  with  local  speed  a,  or  Lagrangian  speed 
0^,  in  directions  of  increasing  and  decreasing  X.  The  characteristic 
parameters  are  chosen  so  that  a>0  in  the  interaction  pulse  and  u<0 
in  the  precursor  simple  wave  region,  while  n=t  on  the  interaction 
front  a  =  0.  In  terms  of  the  variables  q,  y ,  a,  \p  and 


9  t  /  3  t. 

2  =  - —  - — ,  and  the  material  function*  sinH(a) 
3  a/  3  n 

equations  (3.7)  and  (3.8)  are  replaced  by 


=  1  - 


d 

n  A 

d 

n  a 

CS-CD 

2 

C5+CD 


(b 


7^  =  iicose  cH  (2^-  -  cos^  -  qs  i  n  ,  (^.3) 

3a  l  3 n  3n  3  n  J 

7^-  =  Sis  ir.Hcos^  »  (  U  .  It ) 

3  a  3  n 


lL  = 

3  a 


Q  s  i  n  H  s  i  n  t 


-  1  3  a 
1  Dn"  ’ 


(1».5) 


an  d 


2  a 7^-  +  (  2  a-q  cos  ) -^-6  +  sinA-^  =  0, 

3  n  3n  T  3 n 

*  This  relation  is  readily  established  by  using  (  U  .  1  ) 


U.6) 

and  (3.9). 


2  ) 


[2U] 


wh  i  le 

the  definitions  of 

(  a  ■ 

1* 

-  r 

3 1 

3  n 

'  CD 

3  n 

If  tQ 

3t  ,  f  _ 

=  —  and  tn  - 
3  a 

3  t 

—  are 

3  n 

3X_  r  il 

3a  " ' B  3a' 


( 1* .  T ) 


they  are  related  by  the  compatibility  conditions 


3_b.ri  "  iLLa  * 
3  a  3  n 


(  U .  8  ) 


while  elimination  of  X  from  ( ** .  7 )  implies  that 

2aila  +  sinH^ta  +  ffs)=  0. 

In  what  folios  we  take  ( U .  3 )  -  (>*.6),  (>*.8),  (>*.9)  and 


( 1* .  9  ) 


ta  =  otn 


(  »*.  10  ) 


as  a  complete  set  of  equations  for  0 , a , q , Y , t 0  , t n  and  These 

to  be  solved  for  »>0  subject  to  the  condition  that  the  front^O 
borders  an  arbitrary  stretching  wave  vhich  implies  that  on  °  =  0 


tn  -  1  =  X  =  *  =  °, 

while  by  (3. IT)  and  (3-36) 

q  (n)  =  q(0,n),  and  ae(n)  =  a(0,n) 


(  1*.  11  ) 


(  1*.12  ) 


v a. ry  so  that 


q  1 

Ho 


[  1  +  cosH  (  aQ  )  1  a#'  , 


71  ii  2JL 

where  y  <  H  < 


(  1*.  13) 


Note  that  the  only  material  function  which  occurs  in  the 

complete  set  of  equations  (k.3)  -  (^.6)  and  (^.8)  -  (I4.IO'  is 

2  2 

sinH(a).  In  the  important  practical  case  when  2 C  /C  can  be 

U  O 

neglected  compared  with  unity*  we  can  take  H=n  /  2  in  these 
equations  which  then  are  identical  for  all  materials.  Of  course, 
even  in  this  limit,  the  determination  of  X(3,n)  from  (k.7),  and 
the  placing  of  the  characteristics  in  (X,t)  space  involves  the 
material  function  A(a).  However,  at  constant  X  conditions  (  h .  7 ) 
imply  that  n  varies  with  a  so  that 

=  Q(a,n).  (  U  .  1U  ) 

3  a  1  X 

Accordingly,  the  variations  with  t,  at  constant  X,  of  the  variables 
in  (  *4 .  3  )  -  (^.6)  and  (  b  .  8  )  -  (b.10)  can  be  determined  if  one  point 
on  the  (n,a)  curve  defined  by  (U.lt)  is  known;  this  could,  for 
example,  correspond  to  the  arrival  time  n  of  the  front  a=0  at  X. 
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U . 2  Acceleration  front. 

The  expansion  procedure  used  to  determine  conditions  in  the 

pulse  is  best  motivated  by  considering  conditions  at  the  front 

a  =  0  whose  trajectory  is  given  by  (3.2U)  with  a  =  o  and  [ X0 ( 8  )  , 1 0  ( 0  )  ] 

determined  by  conditions  in  the  precursor  wave  by  (3.25)  -  (3. 27). 

At  the  front  the  dependent  variables  in  (  It .  3  )  -  (^.6)  and  ( U . 8 )  - 

( U . 1 0  )  are  continuous,  and  so,  according  to  (  i  .  3  )  -  (  i  .  6  )  ,  are 

-^-5-  and  which  are  given  by  (  1* .  3 )  -  (^.6)  and  (  h  .  1  3  )  as 

3a’  3a  9a 

If  =  notaniH(«o)a;,  §f  =  I!,  s  i  „H  (  a.  )  a;  ,  andff.O  (1.15) 

where,  by  ( *t .  9  )  ,  (^.11)  and  (  ^  .  1  3  )  ,  (  n  )  =  N  0  ,  n  )  =  t  a(  0  ,  n )  varies 

so  that 

!2„exp  [/2 j  r~'  sinjH  (r  ]s  in^  (H(  r )  +  y-)dr]  =  constant,  =  t  say.  (^.16) 

Note  that  flD,  regarded  as  a  function  of  a.,  is  proportional  to  t,  , 
wh  ch  satisfies  (3.26).  41  is  the  only  dependent  variable  which  may 

have  a  discontinuous  first  derivative  at  the  front.  If  (  i  .  6  )  is 
differentiated  with  respect  to  a,  and  the  results  (U.15)  are  used, 
then  an  ordinary  differential  equation  for  |~^-(0.n)  =  ^-(  0  ,  n  )  is 
obtained  which  integrates  to  give 

ao  j 

^•(O.n)exp  [/  r _  cos22H(r)dr]  =  constant,  =  k  say.  (h.17) 

a  Cl 

In  the  special  case  when  H  y-,  (U.l6)  and  (i.17)  simplify  to 
give 

?  3  <f> 

a090  =  c>  and  a— — (0,n)  =  k.  (*t.l8) 

0  3  a 

The  results  (  U  .  1 6  )  -  ( it  .  1 8 )  can  be  used  to  compute  and  compare 

the  values  of  ^-(  t  ,  X  )  and  rf(t,X)  at  the  passage  of  the  front.  These 
3  t  d  L 

are  given  in  terms  of  ac(n)  as 
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ff  =  (2ft)-1  ^(0,n),  and  |f  =j  ( l  +  tan^H(  ae  )  )a^  ,  (*4.19) 

dt  u  d  a  d  L 

which  when  H  t-  can  be  simplified  to  give 

|f  =  ike-1at  and  |f  =  el'0  .  (*4.20) 

Conditions  (*4.19)  and  (*4.20)  imply  that  if  the  stretching  wave 
extends  (contracts)  the  plate  as  it  passes  the  front,  so  that  aQ 
increases  (decreases),  the  I  If  I  increases  (decreases).  In  particular 

a  L 

I  — f-|  -*■  0  as  the  plate  is  unloaded  (a  -*0  )  . 

'  3 1  1  r  o 

However,  since  at  the  front 

1 

|f  =  -X(ao)aQ1  |f  ,  =  -jke  X(ao)aQ  when  H  i  j  ,  (*4.21) 

the  low  frequency  theory  predicts  that  as  a^O  the  curvature  of  the 
plate  become  unbounded.  Again  this  is  prevented  by  the  formation  of 
a  layer,  the  'bending  layer',  in  which  the  low  frequency  membrane 
approximation  is  invalid  and  bending  forces  are  important.  In  this 
layer  the  stretch  varies  across  the  plate  and  as  <*>^0  the  ratio  of 
plate  curvature  in  the  layer  to  plate  curvature  outside  is  0(w  ). 
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U .  3  First  order  solutions. 

Conditions  ( U . 1 9  )  and  ( U . 2  0  )  show  that  at  the  front,  at  times 
when  (a|)  '  is  bounded,  ^7  —  can  be  made  arbitrarily  large  by  taxing 
e /k  sufficiently  small.  Motivated  by  this  result,  and  by  previous 
work  on  finite  amplitude  pulses  by  Varley  and  Cumberbatch  (lU),  (15), 
Varley  and  Rogers  (  l6  )  ,  and  Seymour  and  Varley  (11),  we  seek  asymptotic 
solutions  to  (  U  .  3 )  -  (U.6)  and  (U.8)  -  (U.10),  which  satisfy  (i4.ll)  - 
(U.13)  for  arbitrary  a0(n),  which  are  of  the  form 

N  N  n 

a-ao(n)  =  l  e  a^(a,n),  q-q  (n)  =  [  e  q  ( a ,  n  )  , 

1  °  1 

N  N  N 

y  =  l  e%  (a,n),  t-n  =  l  ent  (a,n),  ft  =  £  €.nft  (a,n)  (U.22) 

1  1  n  1 

N  n 

and  $  =  <t>  (a,n)  +  V  e  <p  (a,n). 
o  L  n 

1 

These  solutions  describe  deformations  in  which  £  varies  slowly 
compared  with  y.  This  is  best  seen  by  noting  that  at  constant  X 


3t /3t  3n  /  3a  3n  ’ 


(U.23) 


which,  together  with  (  U .  3 )  and  (U.22),  implies  that  j^-J  =  0(e) 
as  e-*-0,  (except,  perhaps,  at  isolated  points). 

If  (U.22)  is  inserted  in  (U.6)  this  integrates  to  give 


tan4ii<  =  tan^4  =  A(a  )G(a), 


(U.2U) 


where  the  signal  function  G(a)  is  arbitrary,  and  the  amplitude 
modulation  function 

a° 

A ( a0 )  =  exp[ -/  r  1 c os  2 (  r ) dr ] , 

j  (U.25) 

=  a  ?  when  H  =  t • 

0  c 


Note  that  the  zeroth  order  approximation  ( 1* . 2h  )  gives  the  exact 

variation  of  — ^(0,n)  for  all  a0  ,  when  K=2G'(0).  If  now  {h.22)  is 

3  a 

inserted  in  (  U  .  9  )  «  and  (  U  .  3  )  and  (  U  .  1 U  )  are  used,  an  ordinary 

at. 

differential  equation  is  obtained  for  the  variation  of  9  j  -  ^  a  1  on 
the  a=constant  wavelets:  this ■ integrates  to  give 

a 

o  1  2  2  2  l 

tp  [  -  J  /2r  {s  i  n^H  (  r  )  s  i  n4  (  H  (  r  ) +y)  +  72  cos  2  2  H  (  r  )  .  A  2  (  r  )  G  (a)[l  +  A  (r)G  (a)]  }  dr  ^ 


=  aQ  1  [  1+ a~  1  G2  (  a  )  ]  =  aj'sec2!^ 


when  H 


(h.26) 


terms  of  fl,(n,a),  to  a  first  order. 


(t-n)  =  /n  (n,s)ds, 

O  1 


(14.27) 


1  12 

a  [a  +  a  /g  ( s  )  ds  ]  when  H  = 


while,  according  to  ( b . 7  )  , 


A (aQ) ( X ( 0 , n ) -X )  =  a0(t-n), 


where 


X  ( 0  ,  n  ) 


(I4.28) 


Note  that  n  and  X(0,n)  are  related  to  conditions  in  the  precursor 
wave  on  a=a  =  0  by  (3.2**)  so  that 


n  =  t  o  (  6  )  ,  and  X(0,n)  =  Xo(e). 


(14.29) 


The  statements  (  ^  .  ?  U  )  -  (  U  .  2  8 )  provide  a  first  approximation  for 
the  variation  of  il(X,t)  in  terms  of  the  arbitrary  functions  a^(n)  and 
G(a).  If  (I*. 2U)  is  inserted  in  (3.36)  and  (3.37),  and  if  it  is  noted 
that  to  a  first  approximation  y  :  0 ,  these  then  provide  a  first  approxi¬ 
mation  for  (v,w)  and  (yj,v2).  In  the  special  case  when  ag(n)=0  t!.': 
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first  'approximation'  is  an  exact  solution  and  describes  the  pure 
deflection  waves  of  section  (3.3).  Note,  however,  that  if  a^(n)^rO 
the  expansion  (U.22)  is  only  appropriate  for  a  pulse.  For  at  a 
particle  X  the  assumed  expansion  for  ei,  together  with  (U.lU)  which 
shows  that  x  =  0(e),  implies  that  the  expansion  only  describes 
conditions  over  a  time  interval  after  the  passage  of  the  front  in 
which  a  changes  by  0(e).  The  variation  of  ip  in  this  time  interval 
is,  however,  unrestricted.  Decause  the  expansion  (U.22)  only 
describes  conditions  at  constant  X  over  a  limited  time  it  cannot, 
in  general,  be  used  to  relate  a0(n)  and  G(a)  to  the  driving  con¬ 
ditions  (2.13)  at  X=0.  In  a  future  paper  Parker,  using  techniques 
analogous  to  those  described  in  recent  work  of  Whitham  (IT)  and 
Luke  (9),  will  discuss  high  frequency  expansion  procedures  which 
allow  a  finite  time  variation  at  constant  X. 

r 

If  the  plate  is  only  stretched  so  that  D/_ < 1 ,  which  corre- 

-2  ^s 

spends  to  ,  and  so  that  a.  a#'  is  bounded,  the  expansion  (U.22) 

is  valid  for  all  t_  in  some  region  behind  the  front  a=o ,  and  (U.2U)- 
(U.28)  describe  how  a  bending  pulse  is  modulated  by  the  passage  of 
loading  and  unloading  waves  shich  vary  a# .  Condition  (U.2U)  implies 
that  as  the  plate  is  loaded  (unloaded)  at  an  a=constant  wavelet, 
which  occurs  when  a  increases  (decreases),  |tan24/|  decreases 
(increases)  while,  to  a  first  order,  the  magnitude  of  the  angular 
velocity  of  the  plate  at  a  particle 


3 

3 1 


(2efij 


'  3a 


=  e  '  a0T[  a  ♦G2(u)]  2G'(a)  when  H  =  j-  , 


(  U .  30  ) 


increases  (decreases).  In  particular,  as  the  plate  is  unloaded 

( a.  -»Q  )  ,  (U.2U)  and  (  U .  30  )  predict  that  and  that  -»0  .  However, 

since  in  this  limit  the  filament  theory  predicts  that  and  the 

3  X 

curvature  of  the  plate  are  unbounded  it  is  locally  invalid.  A  bending 
layer  forms  in  which  the  radius  of  curvature  and  thickness  of  the  plate 
are  of  the  same  order  as  o>+0. 


Sc  far  we  have  considered  how  a  deflection  pulse  is  modulated 

by  a  stretching  wave,  and  have  shown  that  to  a  first  order  current 

conditions  at  a  wavelet  a  are  determined  by  conditions  at  a  at  any 

one  previous  time,  and  by  current  conditions  in  the  stretching  wave 

at  the  front  a=0.  To  a  first  approximatio  the  deflection  also 

interacts  with  and  changes  the  mode  of  propagation  of  the  stretchin 

wave.  This  interaction  is  illustrated  by  comparing  the  Lagrangian 
dX 

speed,  —  of  wavelets  carrying  constant  values  of  a  with  the 
Lagrangian  speed  Cg  of  the  B=constant  characteristic  wavelets:  in 
the  precursor  wave,  and  whenever^  =0 ,  these  are  equal.  For  since 


a 


C  (£2- 
da 


a 


-R 


» 


( u. 31 ) 


(**•3),  (1*.13)  and  (U.22)  imply  that  to  a  first  approximation 
CS1frl  a*  l-2sin2^4'.tan^H(l  +  tan^H)_1  (s  in2$*+!taniH[  1-tan^H  ]_1  (  U  .  32  ) 

=  0  when  H  -  y,  UfO), 

=  1  when  H  =  it  . 


Q 

According  to  (U.32),  if  D/_  <1  in  the  bending  pulse  the  wavelets 

carrying  constant  values  of  a  are  slewed  relative  to  the  B=constant 

C 

>1.  When  D/  =1  the  a-constant  wavelets , 

S  CS  “ 

the  B=constant  wavelets  and  the  a=constant  wavelets  move  with  the 

same  Lagrangian  speed  C^. 


wavelets  though  C~  *  77 


—  I 

dt  I 


.  3  Higher  order  approximations 
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Further  terms  in  the  expansions  (U.22)  can  be  obtained  by 
applying  a  simple  iteration  procedure  to  (U. 3)_(U. 5)  and  (U.8)- 
C*4.10).  We  insert  the  expressions  ( U .  3 )  and  (U.10)  for  ^  and  ta 
into  (U.9)  and  regard  this  as  an  ordinary  differential  equation  (I) 
for  fl  on  the  o=constant  wavelets  with  coefficients  depending  on 
a,q  »y  « 4>  »t^  and  their  r,-derivati  ves  .  The  nth  order  approximation 
for  fl  then  satisfies  the  equation  obtained  by  replacing  the  coeffi¬ 
cients  in  (I)  by  their  (n-l)th  approximations.  The  nth  order 
approximations  for  a,q,y  and  t  are  obtained  by  replacing  ft  by  its 
nth  order  approximation,  and  its  coefficient  by  its  (n-l)th  order 
approximation,  and  its  coefficient  by  its  (n-l)th  order  approximation, 
in  the  right  hand  sides  of  (U.3)-(U.5)  and  (U.8):  the  resulting 
equations  are  then  integrated  from  the  front  a=0  along  the  n=constant 
wavelets.  Finally,  the  nth  order  approximation  for  $  satisfies  the 
ordinary  differential  equation  obtained  by  replacing  the  coefficients 
in  ( U . 6 )  by  their  (n-l)^^  approximations.  We  note,  but  do  not  further 
discuss,  that  according  to  ( U . 3 )  this  iteration  scheme,  and  the  formal 

expansions  (U.22),  must  be  modified  if  H-*"  which  occurs  when  C  -*-C  . 

D  b 


As  an  illustration,  which  is  physically  important  and  also  avoids 

IT 

messy  algebra,  we  apply  the  iteration  scheme  when  H=2*  to  calculate  the 
first  approximation  for  the  'slow'  variables  a,q  and  y  in  the  pulse. 
These  approximations,  together  with  the  results  already  established, 
give  the  first  approximations  for  the  time  variations  of  all  dependent 
variables  at  a  particle  X  as  the  pulse  passes.  In  particular,  they 
give  the  first  approximation  to  the  strain  rate  We  obtain, 

to  a  first  approximat i on t 


1  j  (J  j 

a  =  aQ  ♦  aQ  a&  { a  ♦  3ao  /  G  ( s  )  ds  }  , 

O 

q  =  a0  +  ao  1  aj,  |  a  -  a~  /G  (  s  )  ds  }  , 


( U. 33) 


and 


a 

o/G(s ) ds . 

O 


Y=  2ea0 


[  33] 

The  results  (  I*  .  VJ )  indicate  that  our  formal  iteration  scheme 
can  only  be  valid  in  a  puise;  that  is  for  |ea£a|<'"l.  Note  that 
whereas,  to  a  first  approximation,  the  variation  in  the  fast 
variable  ^  at  an  a=constant  wavelet  is  independent  of  conditions 
at  all  precurs’  r  wavelets,  in  any  finite  mp  1  i  t  ud  e  disturbance  the 
slow  variables  are  influenced  by  c  on  a :  t  1  on  ;  at  all  precursor  wave¬ 
lets. 


[34] 


5 .  Shock  J&£I 

In  section  3.1  we  pointed  out  that  according  to  the  low 
frequency  or  long  wave  theory  a  stretching  simple  wave  can  develop 
limit  surfaces,  which  ere  envelopes  of  6-wavelets,  on  which  strain 
gradients  and  strain  rates  become  unbounded.  However,  shortly 
before  this  occurs  the  lateral  velocities  become  comparable  with 
the  longitudinal  velocities  so  that  the  underlying  assumption  of 
that  theory  is  violated.  In  analogy  with  the  laminar  bore  of 
hvdraulics  a  ’shock  layer'  is  generated  in  which  the  strains  and 
velocities  vary  rapidly.  This  layer  is  driven  forward  to  seperate 
two  low  frequency  regions.  By  introducing  a  new  expansion  procedure 
to  describe  regions  of  rapidly  varying  strain,  we  can  obtain  the 
equations  governing  the  structure  of  such  a  layer  and  can,  in 
general,  relate  its  speed  of  propagation  to  the  Jump  in  strain 
ac ros  s  it . 

In  a  layer  where  the  strains  vary  rapidly  (compared  with  their 
values  in  the  low  frequency  region)  the  weighted  scales  (2.1U)  are 
not  appropriate:  the  right  hand  sides  of  (2.16),  (2.17)  and  (2.20) 
are  not  vanishingly  small  as  u-*0 .  If  the  shock  layer  is  centred 
about  the  shock  trajectory  X=S(t)  then  in  terms  of  the  shock  layer 
variables  Y,  t  and 


X  =  u)-1  (X-S(t)  ) 


(5.1) 


equations  (2.16),  (2.17)»  (2.19)  and  (2.20)  become 

a 


ST>2 


-v  aT  3  V 

f^T  .sMt),  )  ♦  : 


a  y 


at 


|T(vi+s*(t)pii)  =  U.-ILL 


(5.2) 


aP 


2  1 


v~+S* ( t )p2 i ) 


(5.3) 
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(5.M 


In  addition  to  (5.2)  -  (5.^),  the  st. esses  are  related  to  the 

deformation  gradients  by  (2.7)  and  satisfy  the  stress  free  boundary 
conditions  =  Tj  =  0  on  Y  =  *1.  As  X  »<•>(-<*>)  the  solutions  to 

( 5 . 2 ) - ( 5  .  b  )  tend  to  the  solutions  of  the  low  frequency  equations  as 
X-»S(t)  from  X>S(t)  (<G(t)),  see  Van  Dyke  (196M). 

Conditions  (  5.3),  together  with  the  boundary  data  at  X  =  +  00 , 
imply  that  to  the  zeroth  order  in  the  shock  layer 

v  +  S'p  =  g  (  t  )  ,  and  v  +  S '  p  =  g  ( t  )  (5.5) 

1  F  1  1  R  1  ’  2  F 2  1  ^2 

so  that  at  any  time  +  S'p^  and  v„  +  S'p2]  are  continuous  across 
the  layer.  To  this  order,  (5.5)  imply  that  the  strain  field  satisfies 
the  compatibility  conditons  (5.M  and,  from  (5.2),  the  equations 


o  .  2 

■S  p 


=  ^r-(T  -S'  2p  ,  )  +  - '■ 

3 X  2  1  F2  1  3  y 


=  0  ,  (5.6) 


which  are  identical  to  the  equilibrium  equations  of  finite  elasticity 

3  p  3p 

with  a  body  force  -S'  ( ^ - ).  If  the  divergence  theorem  is 

applied  to  (5.6),  and  if  the  boundary  conditions  at  the  stress  free 
surfaces  Y  =  ±1  and  those  at  X  =  t®  are  used,  we  obtain  the  result 
that  throughout  the  layer 

y1  [  T]  t (X  ,  r,t )  -  S  ’  2P ,  j  (X  ,r,t ) ) dr  =  f  (t). 


(5.7) 


/  [T21 (X,r,t)  - 


f2  ( t  )  . 


S ' 7  p? )  ( X  ,  r  ,  t )  ]  dr 


(36] 


If  [  v  j  ]  ,  (v  ]  ,  ...  dei.  ote  the  differences  in  values  of  v  j  ,  v . 

at  X  =  ®  and  X  =  -<»,  and  if  it  is  remembered  that  these  differences  are 
independent  of  Y,  (5-5)  and  (5-7)  imply  that 

(  v  j  ]  +  E'[pnl  =  [  v?  ]  +  S’  l  P  p  i  1  =  °>  (5.8) 


and  that 


[ T,  j  ]  +  S'  [ v,  ]  =  [ T2 ,  ]  +  S' l v2  ]  =  0.  (5.9) 

Conditions  (5.8)  and  (5-9)  are  trivially  satisfied,  with  no 
restriction  on  S',  if  (5.6)  have  solutions  which  are  symmetric  about 
X  =  0.  Such  solutions  describe  s  c 1 i t  a  ry  waves .  After  the  passage  of 
a  solitary  wave  conditions  at  a  particle  return  to  those  which  existed 
prior  to  its  arrival.  Note  that  if  the  speed  of  such  a  wave  is  constant, 
solutions  to  (  5  .  U  )  and  (5.6)  satisfy  (  5 . 2  )  -  (  5  .  ^  )  exactly. 

When  the  passage  of  a  shock  layer  cnanges  conditions  at  a 
particle,  (5.8),  (5-9)  and  (3.6)  relate  conditions  in  the  low 

frequency  disturbance  which  occurs  immediately  after  the  passage  of 
the  shock  layer  to  S'  and  conditions  in  the  low  frequency  disturbance 
which  existed  immediately  before  its  arrival.  The  shock  layer  is  a 
narrow  region  in  which  the  low  frequency  approximation  is  locally 
invalid.  The  propagation  speed,  given  in  terms  of  the  Jumps  in 
physical  variables  by  (5.8)  and  (5.9),  is  not  that  of  a  thermodynamic 
shock;  in  the  layer  dissipation  forces  remain  negligible  and  the 
entropy  is  unchanged.* 

If  ve  are  not  interested  in  the  structure  of  the  shock  layer 
but  only  in  h.  it  affects  the  deformation  in  the  low  frequency 
region,  we  may  regard  this  shock  layer  as  a  discontinuity  surface 


*  The  deformation  in  the  shock  layer  may,  however,  produce  velocity 
gradients  which  are  large  enough  for  dissipation  to  be  important. 
Then  thermodynamic  shocks  may  be  formed  or,  in  .he  language  of 
hydraulics,  the  laminar  bore  may  degenerate  into  a  breaking  bore. 


[  37] 


across  which  the  Jump  conditions  (  5  .  8 )  and  (5-9)  hold.  Then 
(5.9)  state  that  the  change  in  force  balances  the  change  in 
momentum,  while  (5.8)  state  that  the  position  of  a  particle  is 
unchanged  by  the  passage  of  the  shock.  These  Jump  conditions  are 
identical  with  those  defined  by  the  conseivation  equations  (3.2) 
and  (3.3)  which  govern  the  low  frequency  deformation.  This  result 
was  anticipated  by  Whitham  (l8)  in  his  work  on  magnetohydrodynamic 
wave’..  He  also  noted  that  the  use  of  the  Jump  conditions  defined 
by  the  conservation  laws  (3.2)  and  (3.3)  can  only  be  Justified  by 
either  having  a  model  for  the  shock  structure,  or  using  sound 
physical  arguments.  Equations  (3.2)  and  (3.3)  can  be  rewritten  in 
many  alternative  ways  as  conservation  laws  and  these  would  define 
Jump  conditions  which,  in  general,  would  be  incorrect.  If,  for 
example,  t  and  the  Eulerian  co-ordinate  x  of  a  particle  are  used  as 
independent  variables  and  the  equations  are  written  in  the  usual 
form  as  conservation  laws  the  associated  Jump  conditions  are 
incorrect.  (In  hydraulics  the  hydraulic  Jump  conditions  are 
derived  from  sound  physical  arguments  and  not  from  the  conser¬ 
vation-type  differential  equations  governing  the  flow,  (see  Stoker 
(12  )  )  . 

There  are  two  possible  ways  of  satisfying  the  Jump  conditions 
(5-8)  and  (5.9).  The  longitudinal  shock,  across  which  ( ip ,  w )  are 
continuous  but  (X,v)  are  discontinuous,  for  which 


S' 


(5.10) 


Such  a  shock  can  be  generated  from  an  initially  shockless  pure 
stretch  simple  wave.  The  transverse  shock,  across  which  X  is 
continuous  but  (v,w,ij»)  are  discontinuous,  which  moves  with 
characteristic  speed 


S  ' 


(5.11) 


[38] 


Th i 3  chock  cannot  be  generated  from  an  initially  shockless  pure 
bending  wave  and  must  Le  regarded  as  the  mathematical  limit  of  a 
non-distorting  simple  wave. 
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1  1  SUPPLEMENTARY  NOTES 

1.'  SPONSORING  Mil  I  '  Ain  activity 
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Disturbances  produced  by  the  motion  of  a  driver  which  is  rigidly 
bonded  to  the  edge  of  a  plate  are  used  to  motivate  parameter  expansion 
techniques  which,  when  applied  to  the  equations  of  finite  elasticity, 
generate  approximating  equations  which  describe  low  frequency  deflection 
and  stretching  waves  travelling  along  stretched  elastic  plates  and  rods 
in  the  limit  when  bending  forces  are  negligible  compared  with  membrane 
forces.  The  structure  of  the  boundary  layer  at  the  driven  edge  and 
shock  layers,  where  the  low  frequency  or  filament  approximations  are  lo¬ 
cally  invalid,  are  also  discussed. 

The  low  frequency  equations  are  used  to  discuss  the  interaction 
between  progressing  finite  amplitude  deflection  and  stretch  waves  in 
the  limit  when  the  stretch  rate  is  small  compared  with  the  angular  speed 
of  the  plate.  The  disturbance  is  locally  that  of  a  pure  deflection 
simple  wave  whose  amplitude  and  frequency  are  modulated  by  slow  varia¬ 
tions  in  the  stretch.  As  the  stretch  increases  the  frequency  increases 
while  the  amplitude  decreases.  The  stretch  wave  is  also  modified  by  de¬ 
flection  of  the  plate:  the  speeds  of  wavelets  carrying  constant  values 
of  stretch  are  always  less  than  their  values  in  the  pure  stretch  simple 
wave  . 
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